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ABSTRACT 

Chandra X-Ray Observatory (CXO) ACIS-S spectra of three quasars which lie behind three foreground 
damped Lya (DLA) absorbers are analyzed in order to attempt to determine the amount of photoelectric 
absorption due to metals present in their x-ray spectra. These absorbers are the three largest neutral 
hydrogen column density absorption-line systems known at low redshift (0.313 < z a bs < 0.524). They 
have HI column densities which lie in the range 3 x 10 21 < Njji < 5 x 10 21 atoms cm~ 2 . At these 
redshifts the amount of photoelectric absorption at x-ray energies is primarily an indicator of the oxygen 
abundance. Since the column densities of these systems are so high, one would expect accurate metailicity 
measurements of them to yield a robust estimate of the column-density-weighted cosmic neutral-gas- 
phase metailicity at z w 0.4. We consider cases where the DLA gas has solar element abundance ratios 
and ones with the a-group element abundance ratios enhanced. For the adopted assumptions, the 
column-density-weighted cosmic neutral-gas-phase metailicity of the non-enhanced elements (e.g. Zn) at 
z w 0.4 likely lies in the range w 0.04 - 0.38 Z Q . 

Subject headings: galaxies: abundances — galaxies: evolution — galaxy formation — nucleosynthesis 
— quasars: absorption lines — quasars: individual (AO 0235+164, S4 0248+430, 
PKS 1127-145) — x-rays: ISM 



1. INTRODUCTION 

Observations of the chemical evolution of the Universe 
provide independent constraints on the processes which 
govern cosmic nucleosynthesis and galaxy formation. This 
includes studies of the chemical makeup of the Universe's 
gaseous and stellar components with increasing redshift. 
With regard to the gaseous structures, it is widely rec- 
ognized that intervening damped Lya (DLA) absorption- 
line systems in quasar spectra are important probes of 
galaxy formation. They are not only the highest neutral 
hydrogen column density systems, with Njji > 2 x 10 20 
atoms cm~ 2 , but consideration of them alone can be used 
to trace nearly all of the neutral gas in the Universe out 
to at least redshift z — 3.5. The work on DLA systems 
has included studies of their cosmological incidence, mass 
density, and column density distribution (e.g., see Rao & 
Turnshek 2000, hereafter RT2000, and references therein), 
the types and impact parameters of the luminous galaxies 
associated with them (e.g., Le Brun et al. 1997, Turnshek 
et al. 2001, Nestor et al. 2002, Turnshek, Rao, & Nestor 
2002, Rao et al. 2003) and their kinematic properties along 
the sight-lines to the background quasars (e.g., Prochaska 
& Wolfe 1999). Of special interest to the work presented 
here, Pettini and collaborators (Pettini et al. 1999 and ref- 
erences therein) have pioneered the method of using DLAs 
to measure the cosmic neutral-gas-phase metailicity of the 
Universe. However, it has not been possible to extend this 



work significantly below z sa 0.6 primarily due to the fact 
that, at these redshifts, the measurements must be done 
from space. This is because there are no suitable unsat- 
urated heavy-element resonance lines in the optical which 
can be used to reliably measure metallicities in high- Nhi 
systems. For example, ZnIIAA2026,2062, which is a good 
tracer of the neutral gas and is relatively unaffected by 
depletion onto dust grains, is only visible in optical spec- 
tra at redshifts z > 0.6; and while CaIIAA3934,3969 and 
NaIAA5891,5897 are visible in optical spectra, they are af- 
fected by complicated depletion processes and possible sat- 
uration effects. Moreover, although many measurements 
at z — can be obtained from observations within the 
Galaxy, these measurements may not be representative of 
the Universe as a whole. Therefore, given that neutral-gas- 
phase metallicities have been measured for DLA systems 
at high redshifts, it would seem that the more appropri- 
ate way to empirically establish the mean metailicity of 
the Universe at low redshift is to measure the asymptotic 
trend (and eventually the spread) in metailicity as one 
probes from high redshifts to low redshifts. 

Here we report attempts to use Chandra X-Ray Ob- 
servatory ( CXO) ACIS-S observations of three quasars to 
measure the amount of photoelectric absorption due to 
metals in their spectra arising from three foreground low- 
redshift DLA systems. The three DLA systems studied 
here have 0.313 < z a b s < 0.524; they are the largest HI gas 
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columns presently known at low redshifts (RT2000), with 
3 x 10 21 < Nhi <5x 10 21 atoms cm -2 as inferred from 
their DLA or 21 cm absorption lines. Since the K shell ion- 
ization edge of oxygen occurs at a rest frame energy of 0.53 
keV and oxygen has relatively high absorptivity in com- 
parison to the other heavy elements (Wilms, Allen, & Mc- 
Cray 2000), the amount of photoelectric absorption in the 
observed 0.3-8 keV energy range, where CXO ACIS-S is 
the most sensitive, is primarily an indicator of a DLA sys- 
tem's oxygen abundance. Furthermore, because the three 
systems are such high-iVjjj systems and the sum of their 
column densities are similar to those probed at higher red- 
shifts, definitive measurements of the metal abundances 
in them would provide a robust determination of the cos- 
mic neutral-gas-phase metallicity at z w 0.4. This is an 
important constraint on the process of cosmic chemical 
evolution; specifically, it would be a measurement of the 
metallicity of the neutral gas which is an important con- 
stituent of matter in the cyclic process of galaxy formation 
(e.g., collapse of gas to form large neutral gas structures, 
followed by the formation of molecular clouds and star 
formation, eventually leading to metal enrichment of the 
environment, and so on — but see Allen 2002 for a different 
perspective). 

In principle, the use of photoelectric absorption to 
measure DLA metallicities has advantages. First, since 
atoms tied up in all states (e.g., including ionized species, 
molecules, and grains) contribute to photoelectric absorp- 
tion, it can provide a more direct measurement of the to- 
tal metallicity. Ionization corrections and depletion onto 
grains need not be considered. Second, due to shielding be- 
low the hydrogen Lyman limit in the absorber rest frame, 
the ionized hydrogen gas fraction is expected to be mini- 
mal relative to the total column, and searches for molec- 
ular hydrogen in DLA systems usually yield null results 
or very low molecular hydrogen fractions (e.g., Levshakov 
et al. 2002, Petitjean et al. 2002, and references therein). 
Thus, the HI column densities inferred from the DLA or 
21 cm absorption lines of these systems are indicative of 
the systems' total column densities. One would therefore 
expect x-ray measurements of the DLA metallicities to be 
straightforward, in the sense that it would not be necessary 
to make corrections that depend on physical conditions. 
However, the main limitation of our CXO ACIS-S obser- 
vations is that they are insufficient to readily resolve the 
oxygen k-shell ionization edge at a particular redshift. In- 
stead, any photoelectric absorption signature is due to the 
cumulative effect of all the metals which may be present in 
objects at all intervening redshifts, and it is not possible 
to definitively determine the abundance of any individual 
element, or even isolate the amount of absorption due to 
metals at any particular redshift. 

Our main aim here is to use CXO ACIS-S observations 
with reasonable assumptions to derive a preliminary mea- 
surement of the mean cosmic neutral-gas-phase metallicity 
of the Universe at redshift z m 0.4. We consider the mea- 
surement preliminary for three reasons. First, since ab- 
sorption edges are not specifically identified, any detected 
absorption could arise from gas at some other redshift. 
However, since the probability of encountering a second 
strong x-ray absorber at z > along any of these sight- 

7 See cxc.harvard.edu/cal/Links/Acis/acis/CaLprods/qeDeg/index.html. 



lines is small, we will simply assume that any detected 
absorption originates at the DLA redshifts. If this were 
not true, our measurements would then represent upper 
limits on the DLA systems' metallicities. Second, it is 
possible that the absorbers we have studied here, while 
sampling the largest HI columns known at low redshift, are 
biased against including more chemically evolved systems 
which contain significant dust, thereby dimming back- 
ground quasars and eliminating them from current sam- 
ples that are being studied. This selection bias, if present, 
could cause us to underestimate the mean cosmic neutral- 
gas-phase metallicity at z « 0.4. Third, CXO ACIS-S cal- 
ibrations may change, especially at energies E < 0.4 keV. 
Nevertheless, we believe that the results reported here are 
valuable because they will, at the very least, provide some 
insight on attacking this problem in the future. 

The organization of this paper is as follows. In §2 we 
present the CXO ACIS-S quasar observations along the 
sight-lines through three low-redshift DLA absorbers, and 
the analyses and assumptions which lead to DLA metal- 
licity determinations. In §3 we combine our results with 
the DLA metallicity measurements of Pettini et al. (1999) 
at z > 0.6; this provides a preliminary estimate of the 
evolution of the column-density-weighted cosmic neutral- 
gas-phase metallicity at z < 3.5. A discussion of the limi- 
tations of these results and their implications is presented 
in §4. 

2. OBSERVATIONS, ANALYSES, AND NEUTRAL-GAS-PHASE 
METALLICITY DETERMINATIONS 

Observations of the three quasars were made with ACIS- 
S CCD #7 on CXO (see cxc.harvard.edu) . The journal of 
observations and other details are given in Table 1. First, 
we emphasize that an essential data processing step is to 
take into account the known decay of the ACIS quan- 
tum efficiency with time. 7 Without this step, the derived 
metallicities for two of the DLA absorbers would have been 
about 0.3 dex higher. The final processing of each observa- 
tion was performed using the XAssist software (Ptak 2001 
and see www.xassist.org), which is a software package that 
assists in x-ray data reprocessing, initial source detection, 
and analysis. The data reductions were performed using 
CIAO 2.2 and the data were reprocessed using CALDB 
2.9. Briefly, XAssist performs the basic data reduction 
steps recommended by the Chandra X-Ray Center ( CXC) 
"threads." We also implemented an optional step which 
removed the 0.5 pixel position randomization. Sources 
were detected using the wavdetect software, which detects 
sources using a wavelet procedure, and the background 
light curves were examined, which permitted removal of 
data taken during times of background flaring. Images of 
detected sources were fitted with elliptical Gaussian mod- 
els in order to determine their spatial extent. The spectra 
were then binned to 20 counts per channel so that the \ 2 
statistic could be used to examine the goodness of each 
spectrum fit. 

The DLA absorption redshifts and HI column densities 
for the three systems are included in Table 1. The HI col- 
umn densities for the AO 0235+164 and PKS 1127-145 
DLAs were derived by fitting Voigt damping profiles to the 
Lya absorption lines in their HST UV spectra (RT2000). 
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Unfortunately, the HST UV spectrum of S4 0248+430 has 
very low signal-to-noise ratio, so the HI column density of 
this system had to be derived from the its 21 cm absorption 
profile (Lane & Briggs 2001) for an assumed spin temper- 
ature of 700 K. Adopting T s = 700 K for this system is 
consistent with the measured spin temperatures of other 
DLA absorbers (Carilli et al. 1996; Lane et al. 1998). 

Absorption of these quasars' x-ray spectra is due to hy- 
drogen, helium, and metals in the Milky Way Galaxy's 
ISM and in other intervening systems. To derive DLA 
metallicities we used the publicly-available XSPEC soft- 
ware package (Arnaud 1996 and see xspec.gsfc.nasa.gov), 
which is often used for analysis of x-ray spectra. We have 
limited our analysis to the 0.4 < E < 8 keV region due 
to calibration uncertainties outside this range, especially 
below 0.4 keV. We proceeded as follows. We assumed that 
a power-law with index L fitted to each quasar's spectrum 
at energies 2 < E < 8 keV could be extrapolated to lower 
energies to provide a reliable indication of the unabsorbed 
quasar x-ray spectrum. This is a reasonable assumption, 
since in no case is the absorbing column large enough to 
affect the continuum at E > 2 keV. Thus, the 2 < E < 8 
keV range fixes both the slope and normalization of the 
x-ray spectrum. Our tests using other quasars with simi- 
lar emission redshifts in the CXO archives, which are not 
in environments where absorption at the emission redshift 
may be prevalent and which have no known DLA absorp- 
tion, show that such extrapolations are generally consis- 
tent with quasar spectra in the 0.4 < E < 2 keV range 
after accounting for Galactic absorption. For the three 
quasars considered here, the best-fit derived Ts over this 
energy range are reported in Table 1. We assumed that 
the most accurate HI column densities of Galactic ISM 
gas along the sight-lines to the three quasars are those de- 
rived from the 21 cm all-sky emission map of Hartmann 
& Burton (1997). These Galactic HI column densities are 
specified in Table 1. However, it should be noted that the 
0.5° x 0.5° angular resolution of the HI map is much larger 
than the pencil beam of a quasar line of sight. Therefore, 
the HI column density derived from a 21 cm emission ob- 
servation may only be approximate since the measured 
value is an average over the size of the beam. We also as- 
sumed that the metallicity of the Galactic gas is given by 
the ISM metallicity model of Wilms et al. (2000). Wilms 
et al. present a good justification for the metallicities they 
adopt. The oxygen abundance in their model is 0.58 Z©, 
while the metallicities of nitrogen and carbon are 0.76 Z© 
and 0.60 Z Q , respectively. The Wilms et al. ISM metal- 
licity model is one of the options incorporated into the 
XSPEC software package. At the same time, in order to 
assess how our final results on DLA metallicities depend 
on Galactic-gas metallicity, we also derived DLA metallic- 
ities under the assumptions of zero and solar Galactic-gas 
metallicity, where the solar metallicity is that specified by 
Feldman (1992). We note that recent investigations gen- 
erally indicate that the Galactic-gas metallicity is demon- 
strably sub-solar, i.e., « 0.65 Z Q (e.g., Sembach et al. 1995, 
Roth & Blades 1995, and Meyer, Jura, & Cardelli 1998). 
For one Galactic sight-line, Meyer et al. (1994) reported 
a Galactic-gas metallicity of 0.4 Z Q . 

For the assumptions given above, and for two differ- 
ent cases of relative DLA metal abundances, we then used 



XSPEC to determine the best-fit overall metallicities for 
the three DLA systems. For one case we assumed that the 
relative DLA metal abundances had solar ratios as speci- 
fied by Feldman (1992), and for the other case we assumed 
that the a-group element relative abundances (specifically 
O, Ne, Mg, Si, S, Ar, and Ca) were enhanced by a factor 
of 2.5 relative to the Feldman (1992) solar metal ratios. 
Consideration of the enhanced a-group element case seems 
appropriate since the x-ray measurements should be most 
affected by the oxygen abundance, whereas at higher red- 
shift the most reliable measurements of the overall DLA 
metallicity come from studies of zinc which tracks iron, 
but show little or no evidence for depletion onto grains. 
The need to consider the possible effects of a-group abun- 
dance enhancements is clearly evident from studies of a 
variety of stellar populations (e.g., see Prochaska et al. 
2000 and references therein). The degree of a-group met- 
als enhancement generally depends on the overall metal- 
licity with [0/Fc]= — 0.4[Fc/H] being close to the general 
trend. For example, [Fe/H]= —1 (i.e., a metallicity of 0.1 
Z Q ) is indicative of an enhancement of a-group metals 
by a factor of 2.5 times solar. Since the higher-redshift 
DLA measurements of Pettini et al. (1999) indicate that 
zinc metallicities are typically < 0.1 Z©, considering a case 
where a-group element abundances are enhanced by a fac- 
tor of 2.5 seems appropriate. 

The derived overall DLA metallicities and errors are re- 
ported in Table 2 for the six combinations of assumptions 
about Galactic-gas metal ratios (three different cases) and 
DLA-gas metal ratios (two different cases) which were out- 
lined above. Also, because the error source was quite dif- 
ferent, the error in photon index was not propagated to 
arrive at the quoted errors; but this is an insignificant com- 
ponent to the error budget. This is especially true given 
the various interpretive assumptions we have made (i.e., 
the six combinations of assumptions mentioned above plus 
assumptions put forth in §1). Uncertainties in the Galac- 
tic and DLA HI column densities have also not been pro- 
pogated, but again, these contributions to the error budget 
are small relative to the interpretive assumptions. The \ 2 
statistic per number of degrees of freedom is listed in Table 
2 for each fit. This fit refers to the entire 0.4 < E < 8 keV 
range. These statistics show all of the fits to be fairly rea- 
sonable, so this is an indication that x-ray observations at 
such low resolution are unable to distinguish between the 
six cases. Because of the overwhelming evidence that the 
Galactic-gas ISM metallicity is likely to be close to that 
specified by Wilms et al. (2000), we adopt these results in 
the remainder of the paper. 

The left and right panels of Figures 1 — 3 show the ef- 
fect of DLA-gas absorption on the quasar x-ray spectrum. 
Both panels include a power-law fit in the 2 — 8 keV range 
that is extrapolated to 0.4 keV to obtain the unabsorbed 
quasar x-ray spectrum and absorption by Galactic gas ac- 
cording to the Wilms et al. (2000) ISM metallicity model. 
The right panel shows the best-fit model with DLA-gas 
absorption also included. The results indicate that, for 
the AO 0235+164 and S4 0248+430 absorbers, the overall 
DLA metallicities may be substantial, i.e., > 0.1 Z©, and 
metallicities only a factor of w 2 less than Galactic ISM 
values cannot be ruled out. In §4.1.2 we do, however, dis- 
cuss the possibility that another significant absorber may 
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be present in the S4 0248+430 spectrum. As for the PKS 
1127—145 DLA absorber, for our adopted assumptions we 
find that there is no evidence for absorption due to metals 
at the DLA redshift. In fact, our analysis under-predicts 
the x-ray flux in the 0.4 < E < 2 keV region even when the 
DLA gas metallicity is zero. This can be seen by inspect- 
ing the right panel of Figure 3, where the x-ray spectrum 
is seen to lie slightly above the model fit which has zero 
DLA metallicity. This may indicate a problem with our as- 
sumptions in the PKS 1127—145 analysis, which is further 
discussed in §4.1.3. 

We note that Bechtold et al. (2001) also recently per- 
formed an XSPEC analysis to study the metallicity of the 
DLA absorber in PKS 1127—145 (see both our Figure 3 
and Bechtold et al.'s Figure 3). However, pile up, which is 
a problem with the data (G. Chartas, private communica- 
tion), was not considered in their analysis. We attempted 
to model and remove this effect. While the effects of pile 
up can be mitigated by going off axis, one would have 
needed to be > 5 arcmin off axis to begin to obtain a 
significant reduction in pile up (G. Chartas, private com- 
munication) . 

In the Bechtold et al. analysis, some emphasis was 
placed on exploring the parameter spaces of x-ray spec- 
trum fits, including the possibility of absorption at other 
rcdshifts. Ignoring the fact that pile up was not considered 
in their analysis, their work indicated that, indeed, the 
most plausible redshift for x-ray absorption was the DLA 
redshift (z a b s = 0.313). For Galactic absorption along 
the PKS 1127—145 sight-line they used the observations 
of Murphy et al. (1996) and adopted an HI column den- 
sity of 3.8 x 10 20 atoms cm , which is about 16% larger 
than what we used based on the observations of Hartmann 
& Burton (1997). The metallicity they adopted for the 
Galactic ISM gas was presumably solar. However, when 
making fits to any DLA x-ray absorption, they chose to 
fix the DLA gas metallicity to be either zero (primordial) 
or solar (processed), and then they determined best-fit 
Nhi values for absorption of the x-ray spectrum. They 
found that solar metallicity gas with Nhi = 1.2 x 10 21 
atoms cm~ 2 provided the best fit to the DLA-absorbcd 
x-ray spectrum. Since this is 23% of the HI column in- 
dicated by the fit to the DLA line in the HST UV spec- 
trum (RT2000), they inferred a metallicity of 0.23 Z Q for 
the PKS 1127-145 DLA absorber. However, they also 
showed that zero metallicity DLA gas at z a b s — 0.313 with 
Nhi = 3.3 x 10 21 atoms cm" 2 can adequately absorb and 
explain the x-ray spectrum. But in this case the required 
Nhi value is still only 63% of the HI column derived from 
the HST spectrum. This conclusion is implausible since 
(as Bechtold et al. also note) metal absorption lines are 
observed in the quasar spectrum at the DLA redshift and 
the HI column density is greater than Nhi = 3.3 x 10 21 
atoms cm~ 2 . It would then follow that DLA gas with 
the known column density, i.e., Nhi = 5.1 x 10 21 atoms 
cm~ 2 , cannot have a metallicity of 23% solar and still be 
consistent with the x-ray spectrum. 

The appropriate way to proceed, which is the approach 
we took in our analysis, is to run the XSPEC software by 
fixing the DLA HI column density (since it is known) while 
allowing the overall DLA absorber metallicity to be a free 
parameter. This takes advantage of the fact that fitting 



a Voigt profile to the UV DLA absorption line provides 
a metallicity-independent method of determining the HI 
column density. 

3. COLUMN-DENSITY-WEIGHTED COSMIC 
NEUTRAL-GAS-PHASE METALLICITY 

The integrated HI column densities of 42 high-redshift 
DLA systems in five arbitrary redshift bins (e.g., see Pet- 
tini et al. 1999) is comparable to the integrated HI col- 
umn density of the three systems we have studied here at 
z w 0.4. This is illustrated in Figure 4. The distinguishing 
difference of the summed HI column densities illustrated 
in Figure 4 is that the lowest redshift bin (this work) per- 
tains to metallicity measurements at x-ray energies, while 
the higher redshift bins pertain to metallicity measure- 
ments made at UV/optical wavelengths. From this fig- 
ure it is clear that if metallicity determinations of suffi- 
cient accuracy existed for the three individual low-redshift 
DLA systems studied here, a very interesting result could 
be inferred on the mean eolumn-density-weighted cosmic 
neutral-gas-phase metallicity at z « 0.4. For this work, 
we adopt a definition analogous to the one used by Pet- 
tini et al. (1999) to determine the mean column-density- 
weighted metallicity of Zn. In the absence of an effect 
which would bias our sample (see §1), such a measure- 
ment would represent the actual cosmic neutral-gas-phase 
metallicity. Moreover, such a determination would be ro- 
bust in the sense that future metallicity measurements of 
many lower- Nhi systems would only slightly modify the 
result because of their low weight. 

For the assumption of Galactic ISM gas with the metal- 
licity specified by Wilms et al. (2001) we have determined 
the following from the XSPEC fits. If DLA gas has so- 
lar ratio metals, the column-density-weighted neutral gas- 
phase metallicity at z w 0.4 is 0.19 ± 0.11 Z e . Alter- 
natively, if the DLA gas has a-group metals which are 
enhanced by a factor of 2.5 over solar ratios, we find a 
metallicity of 0.10 ± 0.06 Z Q for the non-enhanced ele- 
ments (e.g. Zn). In addition, if one were to take the view 
that the PKS 1127-145 DLA metallicity results should 
be excluded from this analysis (see §2 and §4.1.3), the re- 
sults become 0.31 ± 0.07 Z Q (for solar ratio metals) and 
0.15 ± 0.04 Zq (for non-enhanced metals mixed with en- 
hanced a-group metals). 

As discussed earlier, depletion onto grains is not a con- 
cern for x-ray derived metallicitics. Thus, it is possible to 
make a comparison between the x-ray derived low-redshift 
DLA metallicity results and those derived from the Znll 
lines in high-redshift DLAs. This is because Zn is not 
significantly depleted onto grains and the columns of Zn 
not in the singly-ionized state are expected to be relatively 
small, especially for the large HI column densities involved 
here (e.g. Vladilo et al. 2001). Since we quote metallic- 
ities for the non-enhanced elements in both cases in the 
preceding paragraph, the results can be directly compared 
to the DLA [Zn/H] metallicities considered by Pettini et 
al. (1999) at high redshift. These comparisons are made 
in Figure 5 (all three DLA absorbers) and Figure 6 (ex- 
cluding the PKS 1127-145 DLA absorber). The effect of 
using a Galactic-gas metallicity other than the ISM values 
of Wilms et al. (2001) can be judged from the range of 
results presented in Table 2. 
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Figures 5 and 6 show that the x-ray observations cannot 
be used to rule out a significant increase in the column- 
density-weighted neutral gas phase metallicity with de- 
creasing redshift by z « 0.4. However, given the vari- 
ous assumptions we have made, and the formal errors in 
our metallicity determinations for these assumptions, the 
results should also not be interpreted as evidence for an 
increase in metallicity with decreasing redshift. 

4. DISCUSSION 

4.1. Comments on Individual DLA Systems 

Of course, analysis of the Chandra data is not limited 
to results on DLA absorption in these quasars. Additional 
analysis on the fields containing these quasars can be found 
in the XA ssist output (see www.xassist.org) . With regard 
to individual DLA systems, various details should be noted 
here for each of the objects studied. These are discussed 
below. 

4.1.1. AO 0235+164 (z abs =0.524) 

AO 0235+164 is a well-known object that has prop- 
erties of so-called optically violent variables (OVVs). It 
has been described as a blazar, a BL Lac object, and a 
high-polarization QSO (e.g., see Burbidge et al. 1996 and 
references therein). The high column density absorption 
system at z a b s = 0.524 was first discovered in redshifted 
21 cm absorption by Roberts et al. (1976). Subsequently, 
Wolfe et al. (1982) showed the 21 cm absorption to be 
variable and Briggs (1983) proposed a kinematic model 
for the absorption. 

Cohen et al. (1999) used their HST UV spectrum 
of AO 0235+164 to study the DLA line associated with 
the 21 cm absorption and deduce the HI column den- 
sity and other properties of the absorber. They found 
Nhi = 5x 10 21 atoms cm' 2 ; they also reported evidence 
for dust at the absorption redshift, but at a significantly 
lower level than what would be expected for a Galactic 
line-of-sight at similar HI column density. Our indepen- 
dent analysis of the HST UV spectrum of AO 0235+164 
yielded Nhi = 4.5 x 10 21 atoms cm~ 2 for the DLA ab- 
sorber, which is the value adopted in our analysis (Tables 
1 and 2). 

Using data taken with the Einstein Observatory Imag- 
ing Proportional Counter in August 1980, Madejski (1994) 
performed an analysis which showed the metallicity of the 
DLA absorber to be 2 + 1 Z for the assumptions he made. 
This is a higher metallicity than what we infer from the 
superior Chandra data. 

Finally, HST imaging and follow-up spectroscopy of the 
field by Burbidge et al. (1996) has led to a rather sur- 
prising serendipitous discovery about the environment sur- 
rounding the AO 0235+164 DLA absorber. There are two 
optical identifications of objects (called "Al" and "A") 
at redshift z w 0.524 that lie within 2 arcsec of the AO 
0235+164 sight-line. The nearer object, Al, is separated 
by 1.08 arcsec from the sight-line and has a WFPC2 in- 
strumental magnitude of mp^ 2W = 21.3. It is a resolved 
emission-line galaxy, appearing as a structure which is 
elongated 2 arcsec, with a brighter inner region that has an 
extent of 0.5 arcsec. Object A is separated by 1.96 arcsec 
from the sight-line and has mF702W = 20.6. It appears as 
an unresolved point source surrounded by faint nebulosity. 



Its spectrum exhibits broad emission lines as well as broad 
absorption lines (BALs), which have a maximum outflow 
velocity (6000 km s _1 ); this maximum velocity is on the 
low side of what is normally observed in BAL QSOs. How- 
ever, Burbidge et al. (1996) note that, based on object A's 
magnitude and the surrounding faint nebulosity, it might 
best be classified as a Seyfert 1 AGN, in which case the 
outflow velocities are rather large in comparison to those 
sometimes seen in AGN. In any case, the presence of both 
an emission- line galaxy (object Al) and a BAL QSO/ AGN 
(object A) in the vicinity of the AO 0235+164 DLA ab- 
sorber is suggestive of a chemically-evolved environment. 

4.1.2. S4 0248+430 (z abs = 0.394J 

The identification of the DLA absorber in S4 0248+430 
is based on 21 cm absorption (Lane & Briggs 2001). As 
such, a DLA absorption line was not used to derive its 
HI column density. Instead, a reasonable value of the 
spin temperature was assumed for the HI column density 
derivation (§2). The fact that the continuum was too faint 
to observe the DLA line in S4 0248+430's HST UV spec- 
trum (RT2000) may be an indication that significant dust 
is present along the sight-line. 

This sight-line also represents a case where significant 
metal-line absorption from objects other than the DLA 
system are clearly present. However, it is unclear if the 
column densities of these other absorption systems are 
sufficient to give rise to x-ray absorption. In particular, 
in addition to the DLA system at z a b s = 0.394, Womblc 
et al. (1990) identify two other distinct low-ionization sys- 
tems at z a bs — 0.052 and z a b s = 0.452. The rest equivalent 
widths of the Mgll lines at z a b s = 0.452 (Womblc et al. 
1990) arc small enough that this system is not likely to 
be in the DLA regime (RT2000). Indeed, our HST spec- 
trum of the Lya line at z abs = 0.452 (RT2000) indicates 
that this system is not damped. Thus, we conclude that 
the 

z abs — 0.452 does not give rise to any appreciable x-ray 
absorption. However, as discussed below, x-ray absorption 
from the z a b s = 0.052 is a possibility. 

The system at z a b s = 0.052 is likely associated with 
a luminous spiral galaxy at redshift z = 0.052 which 
is offset 14.7 arcsec from the sight-line to S4 0248+430 
(Womble et al. 1990; Sargent and Steidel 1990). At such 
low redshift, light from this extended luminous galaxy 
overlaps and contaminates the image of the background 
quasar. Both Womble et al. (1990) and Sargent & Steidel 
(1990) detected Call and Nal absorption at this redshift 
in the quasar's spectrum. The high-resolution spectrum of 
Womble et al. shows two primary absorption components 
at this redshift. Their curvc-of-growth analysis indicates 
that the total Ca+ column density of both components is 
w 5.3 x 10 13 atoms cm~ 2 . If there were no calcium de- 
pletion onto grains and this were indicative of a neutral 
gas column with solar abundances, the total HI column 
density would be « 2.3 x 10 19 atoms cm -2 , which would 
not be large enough to cause significant x-ray absorption. 
However, the curve-of-growth error analysis of Womble et 
al. admits the possibility that the Ca + column density is 
significantly higher. Moreover, calcium is easily depleted 
onto grains and, as discussed above, there is evidence for 
dust reddening along this sight-line. Thus, we cannot rule 
out the possibility that some of the x-ray absorption seen 
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in the S4 0248+430 spectrum comes from this system. 

4.1.3. PKS 1127-145 (z abs = 0.313; 

The properties of the various emission- line galaxies near 
redshift z — 0.313 along the sight-line towards PKS 
1127-145 are described by Lane et al. (1998), Nestor et 
al. (2002), and Rao et al. (2003). Most notably, there is 
a large patchy/irregular low surface brightness structure 
which is very extended in the north-south direction and 
~ 3.5 arcsec to the west of the sight-line. 

As discussed in §2, Bechtold et al. (2001) originally 
studied the possibility of x-ray absorption from this DLA 
system. This same group has also made a detailed study 
of the Chandra x-ray images of the field (Sicmiginowska et 
al. 2002). Recall that the main problem we encountered in 
our analysis of x-ray absorption from this DLA absorber 
(§2) was that, even with zero metallicity DLA gas, after 
Galactic absorption was taken into account the DLA ab- 
sorption column required to explain the x-ray spectrum at 
energies below 2 keV was less than what was measured in 
the RT2000 HST UV spectrum. This suggests that extrap- 
olation of the power-law fit to the 2 — 8 keV region to lower 
energies under-predicts the x-ray flux in the 0.4 — 2 keV 
region, indicating that there may be, for example, a signif- 
icant soft x-ray energy excess in PKS 1127—145. Siemigi- 
nowska et al. (2002) present some evidence for this; how- 
ever, the flux from an x-ray jet identified by Sicmiginowska 
et al. in the x-ray image appears to be negligible. In any 
case, estimating the metallicity of this object is problem- 
atic. We can only conclude that in this DLA system there 
is presently no evidence for relatively high metallicity, un- 
like the evidence in the other two DLA absorbers. 

4.2. Limitations of the Results and Implications 

The details of the objects studied (§4.1) and the spec- 
tral resolution of the CXO ACIS-S observations present 
limitations to the direct detection of metals at the DLA 
rcdshifts. Thus, our results and interpretation depend on 
the reliability of the assumptions we have made. For now, 



we believe that our assumptions represent the most prob- 
able case. This leads to the conclusion that we can not 
rule out the possibility of significant evolution of the cos- 
mic neutral-gas-phase metallicity in the form of increasing 
metallicity with decreasing redshift. In fact, Figure 6 sug- 
gests this trend, but the results are also consistent with no 
evolution (Figure 5). 

The statistics of low-redshift DLAs (RT2000) support 
the assumption that the identified DLA system in each 
quasar spectrum is the dominant large column density sys- 
tem along the sight-line to the quasar. Except in the case 
of S4 0248+430, the possibility of other significant absorp- 
tion can be ruled out with some degree of certainty (§4.1). 
Also, because our sample is a radio-selected one, it proba- 
bly does not suffer from the bias mentioned in §1, namely 
that chemically evolved DLA systems with significant dust 
are excluded from the sample due to optical dimming of 
background quasars. If present, this selection bias would 
cause an underestimate of the mean cosmic neutral-gas- 
phase metallicity at z « 0.4. 

The results reported here show that CXO ACIS-S ob- 
servations can be used to investigate some important is- 
sues involving the x-ray absorption properties of DLA gas. 
However, our present results on the objects considered here 
indicate that x-ray grating observations are required in or- 
der to definitively make DLA metallicity measurements at 
x-ray energies. Such observations would not only permit 
individual features from redshifted metals to be observed, 
it would also eliminate the simplifying assumption that the 
unabsorbed x-ray spectrum can be modeled with a simple 
power law. 
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Table 1 

Journal of Chandra X-Ray Observatory ACIS-S Observations 





(10 21 cm- 2 ) 1 


DLA 

Zabs 


Ngf A Redshifted 
(10 21 cm- 2 ) DLA Oxygen 
K edge (keV) 


Obs. 
Date 


Cleaned 
Exp Time 
(ksec) 


Total 2-8 keV Program 
Counts T 2 ID (PI) 


AO 0235+164 0.94 
S4 0248+430 1.31 
PKS 1127-145 1.19 


0.87 
0.89 
0.33 


0.524 
0.394 
0.313 


4.5+0.4 0.35 
3.9+0.6 3 0.38 
5.1+0.9 0.40 


20/09/00 
27/09/00 
28/05/00 


27.5 
22.8 
27.3 


3,577 1.61+0.10 884 (Turnshek) 
2,570 1.60+0.12 885 (Turnshek) 
15,818 1.20+0.05 866 (Bechtold) 



Note. — (1) Galactic HI column density derived from the 21 cm emission map of Hartmann & Burton (1997). Uncertainties 
are on the order of 2%. (2) The quoted error is approximately a lcr error derived from the x 2 statistics of the fit. (3) The HI 
column density is derived using the 21 cm data reported in Lane & Briggs (2001) and is slightly revised from the Lane (private 
communication) result reported in RT2000. 



Table 2 

Derived Neutral-Gas-Phase Metal Abundances for Three Low-z DLAs 1 







Galactic 




nSt a 

HI 


Best-Fit DLA 


Best-Fit DLA 


X 2 per dof 4 


X 2 per dof 4 


Vc/ LiCL kJ CXii. 


(10 21 atoms Metallicity 




( 1 21 atoms 


Mctallicitv (ZrO 


Metallicitv (7,^) 


f solar ratio 


(prinaripprl rv 




per cm 2 ) 


Model 2 




pel Clll J 




( r\ n Vi o n r- 1~\ r\ vn/af olc^ 






AO 0235+164 


0.87 





0.524 


4.5+0.4 


0.61+0.07 


— U.Uo 


150/131 


149/131 






ISM 






0.37+0.06 


0.19+0.03 


143/131 


143/131 






Solar 






0.24+0.06 


0.11+0.03 


142/131 


142/131 


S4 0248+430 


0.89 





0.394 


3.9+0.6 


0.46+0.08 


0.22+0.04 


86/100 


86/100 






ISM 






0.23+0.07 


0.11+0.04 


84/100 


84/100 






Solar 






10+ 008 

u - iu -0.07 


04+ 04 


83/100 


83/100 


PKS 1127-145 5 


0.33 





0.313 


5.1+0.9 


(< 0.0035) 


(< 0.0024) 


389/300 


389/300 






ISM 






(< 0.0032) 


(< 0.0022) 


367/300 


367/300 






Solar 






(< 0.0029) 


(< 0.0020) 


475/300 


371/300 



Note. — (1) The derived DLA metallicities hold for the specified values of Galactic (z = 0) Nri and metallicity. If additional 
absorption is present, these results represent upper limits on DLA metallicity (see text for assumptions and discussion). (2) Three 
Galactic-gas metallicity models were considered: zero metallicity, the ISM metallicity model of Wilms et al. (2001), and the solar 
metallicity model of Feldman (1992). The ISM metallicity model used to construct Figure 5 has an oxygen abundance of 0.58 Z , a 
nitrogen abundance of 0.76 Z , and a carbon abundance of 0.60 Z Q (see §3). (3) The DLA metallicities quoted here were inferred 
from the XSPEC results so that they could be compared to the zinc metallicity results at high redshift (e.g., Pettini et al. 1999), 
where the DLA zinc abundance is given relative to the solar zinc abundance. In the case of DLA gas with a-group elements enhanced 
by a factor of 2.5, one sees that the inferred equivalent zinc metallicity is lower; this is because the abundance of elements like oxygen 
are enhanced in the gas, allowing the abundance of elements like zinc to be lower and yet produce the best fit to the x-ray spectrum. 
The quoted errors are approximately la for the assumptions made. (4) The \ 2 P er degree of freedom for the fit to the x-ray spectrum. 
(5) Spectrum corrected for pile up. However, the inferred low metallicity may be an artifact of assuming that the unabsorbed x-ray 
spectrum can be described by a single power law. 
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0235+164 0235+164 




Energy (keV) Energy (keV) 

Fig. 1. — Chandra ACIS-S spectrum of the quasar AO 0235+164. The left panel shows the power-law fit to the 2 — 8 keV region of 
the spectrum, extrapolated to 0.4 keV while only accounting for Galactic ISM absorption. The right panel shows the same power-law fit 
extrapolated to 0.4 kev but including the effects of both Galactic ISM absorption and DLA absorption. Details are reported in the text and 
tables. 
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0.5 




Energy (keV) 



Energy (keV) 



Fig. 2. — Chandra ACIS-S spectrum of the quasar S4 0248+430. The left panel shows the power-law fit to the 2 — 8 keV region of 
the spectrum, extrapolated to 0.4 keV while only accounting for Galactic ISM absorption. The right panel shows the same power-law fit 
extrapolated to 0.4 kev but including the effects of both Galactic ISM absorption and DLA absorption. Details are reported in the text and 
tables. 
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1127-145 1127-145 
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Energy (keV) Energy (keV) 

Fig. 3. — Chandra ACIS-S spectrum of the quasar PKS 1127—145. The left panel shows the power-law fit to the 2 — 8 keV region 
of the spectrum, extrapolated to 0.4 keV while only accounting for Galactic ISM absorption. The right panel shows the same power-law 
fit extrapolated to 0.4 kev but including the effects of both Galactic ISM absorption and DLA absorption. Note that fixing Nhi in the 
DLA system to the value derived from the fit to the DLA line in the HST UV spectrum and assuming zero DLA metallicity causes us to 
under-predict the flux in the x-ray spectrum in the 0.4 < E < 1 keV range. Details are reported in the text and tables. 
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Fig. 4. — Individual HI column densities at the observed rcdshifts and summed HI column densities in various redshift intervals. The filled 
star symbols represent the three low-redshift DLA systems studied here; the filled square symbols represent the individual highcr-redshift 
DLA systems studied by Pettini ct al. (1999), Prochaska & Wolfe (1999), Molaro ct al. (2000), and Ge, Bechtold, & Kulkarni (2001). Note 
that some of the data used in Pettini ct al. (1999) were below the DLA threshold column density (not shown here), but these points have 
little affect on the column-density- weighted mctallicity results. The unfilled symbols are the summed HI column densities in the designated 
horizontal redshift bins. 
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Fig. 5. — Log of the column-density-weighted cosmic neutral-gas-phase metallicity derived from DLA systems versus cosmic time (bottom 
horizontal axis) and redshift (top horizontal axis) for a cosmology with Ho = 65 km s _1 , qo = 0.5, and A = 0. The low-redshift metallicity 
determination was derived from XSPEC analysis of the Chandra ACIS-S x-ray spectra of three quasars containing DLA systems (Table 1 
and 2). The horizontal "error bars" represent the redshift intervals covered (e.g., 0.313 < z a i s < 0.524 for the low-rcdshift point). The 
high-rcdshift metallicity determinations shown as open squares are those for zinc (a non-a-group element) as reported by Pettini et al. (1999). 
The two low-redshift metallicity determinations shown here (which are slightly displaced in redshift so that they can be distinguished in the 
figure) are those derived under the assumption that Galactic (z = 0) ISM gas absorption follows the metallicity model recommended by Wilms 
et al. (2001). Other assumptions are discussed in §2 of the text. The point plotted with an open triangle was derived under the assumption 
that the DLA gas has solar ratio metals as specified by Feldman (1992). The low-redshift point plotted with an open circle was derived under 
the assumption that DLA a-group metals where enhanced by a factor of 2.5 over solar ratio metals. By considering these two cases (a-group 
elements not enhanced and a-group elements enhanced), a fair comparison can be made with the higher redshift metallicity determinations 
based on zinc. The horizontal dashed line indicates solar metallicity; the horizontal dotted line indicates Galactic ISM gas metallicity. 
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Fig. 6. — Same as Figure 5, but with the zero metallicity results on the PKS 1127—145 DLA absorber removed. The possible justification 
for removing this DLA absorber from the analysis is discussed in §2 and §4.1.3. 



